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ABSTRACT 

We probe the spatial and dynamical structure of the old open cluster M67 
using photometric data from the Sloan Digital Sky Survey's sixth data release. 
Making use of an optimal contrast, or matched filter, algorithm, we map the 
distribution of high probability members of M67. We find an extended and 
elongated halo of likely members to a radius of nearly 60'. Our measured core 
radius of R cor e = 8'24 ± 0'60 is somewhat larger than that of previous estimates. 
We attribute the larger core radius measurement to the SDSS probing lower mass 
main sequence stars than has been done before for similar studies of M67, and the 
exclusion of post main sequence M67 members in the SDSS sample. We estimate 
the number of M67 members in our SDSS sample to be 1385 ± 67 stars. A lower 
limit on the binary fraction in M67 is measured to be 45%. A higher fraction of 
binary stars is measured in the core as compared to the halo, and the luminosity 
function of the core is found to be more depleted of low-mass stars. Thus the 
halo is consistent with mass segregation within the cluster. The galactic orbit of 
M67 is calculated from recent proper motion and radial velocity determinations. 
The elongated halo is roughly aligned to the proper motion of the cluster. This 
appears to be a result of mass segregation due to the galactic tidal field. Our 
algorithm is run on 2MASS photometry to directly compare to previous studies 
of M67. Decreasing core radii are found for stars with greater masses. We test the 
accuracy of our algorithm using 1000 artificial cluster Monte Carlo simulations. 
It is found that the matched filter technique is suitable for recovering low-density 
spatial structures, as well as measuring the binary fraction of the cluster. 

Subject headings: open clusters and associations: individual (M67) - galaxy: 
kinematics and dynamics 
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Introduction 



For over 40 years the spatial str ucture of s tellar clusters has been characterized 
by the widely accepted King model 



range of observed clusters well (e.g. 



(King 



1962) . This density model matc hes a wide 



McLaughlin fc van der Marel 



20051 ). As 



Kind fl!966h 



notes, however, the model does not describe the specific dynamical evolution or state 
of a particular cluster. As clusters age their constituent stars inevitably undergo close 
gravitational encounters with other members of the cluster. Through energy equipartition, 
lower mass stars are given higher velocities, and thus larger orbits in the cluster. Many will 
be given velocities greater than the escape speed for the cluster's gravitational well. We 
observe this dynamical equipartition in the segregation of masses radially across a cluster, 
with lower mass stars being found preferentially further from the center than higher mass 
stars. Binary stars are also more centrally concentrated than their single-body counterparts. 
The King model does not deal with the unique history a cluster may have within its parent 
galaxy, or the state and location of stars previously associated with it. 

Escaped stars may be projected across large areas on the sky. Within disrupting 
clusters (such as most open clusters for example) the ejected stars may be a significant 
fraction of the cluster's initial mass. It is also likely that the ejected stars are not 
distributed in a spherical manner, due to tidal disruptions from the galactic potential or 
close encounters with giant molecular clouds. By understanding the total amount and rate 
of mass lost in a cluster, we may begin to paint a picture of the cluster initial mass function, 
and the history of its disruption in the galaxy. Since open clusters are plentiful and found 
with a wide variety of ages, we may hope to study this process at many stages. 

M67 (NGC 2682) is an old (~ 4Gyr), nearby (~900pc), and well-studied open cluster. 
Its proximity has made it an i deal target for a wide array of studies over the past centu ry, 



ranging from radial velocities (iMathieu et al. 



19861 ) and proper motions (1Sanderslll977l : 
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Girard et al.lll989f ). to late-type stellar evolution and binary fractions ( [Montgomery et al. 
19931 ). 



The disruption scenari os and timesca 

~3~ 



studied for some time (e.g. 



Spitzer 



es of open clusters in our galaxy have been 



1958ajjb|). Recently N-body models have become 



sophisticated enough to place constraints on th e role of giant molecular clouds and deta iled 



initia l structure in cluster dissolution (e.g. see 



Gieles et al 



2006 



Gieles fc Baumgardt 



2008 1). Since most open clu sters are thought to generally dissipate over ~ 10 7 year timescales 



( iBinney fc Tremaine 



19871 ). the existence of old open clusters is som ewhat of an anomaly 



While M67 is a classic example of a highly evolved open cluster (e. 



g. 



van den Bergh 



Carraro et al. 



1953), 



20061 ) 



older examples do exist. NGC 6791 has an age of at least 8Gyr ((( 
roughly twice that of M67. This open cluster is believed to be one of the most massive 
in our galaxy, and its high stellar density is surely a prerequisite for its survival to such 
exceptional age. Clearly this is an abnormal open cluste r, and its origins un der the typical 



open cluster formation scenarios have been questioned (ICarraro et al. 



20061 ). Despite its 



uncertain origins, NGC 6791 like other older clusters displays a dy namically evolved mass 



function, indicative of mass segregation ( jKaluzny fc Rucinski 



19951 ). 



Due to the paucity of old open clusters such at M67 and NGC 6791, and the large 
number of younger clusters, it is reasonable to assume that a great many open clusters 
existed in the past and have since been destroyed. These lower mass clusters, like the young 
clusters found today, likely had a wide range of masses and numbers of members. Indeed it 
is believed that most stars which are now part of the galactic disk originated in clusters of 



various sizes. These clusters must have had densities that were sufficient 



them to be completely dissociated within about 600 Myr (IBergond et al. 



y low as to allow 



20011 ). 



Ongoi ng work to fu ll y mo del the evolution of M67 over its entire lifespan has been 



promising. 



Hurley et al. 



( 120051 ) have created an N-body simulation of M67 whose result 
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over 4 Gyr of evolution shows reasonably good agreement with observations of the spatial 
distribution and stellar composition of the cluster. These types of simulations allow 
observers to study the cluster's initial conditions and their effect on the present state. M67 
is found in this model to have lost at least 75% of its stellar mass to dynamical and stellar 
evolution, and it is suggested that the cluster will continue to dissipate. 



Mass segregat ion has previous 
properties of M67 



y been suggested as the cause for the observed structural 



Bonatto fc Bical (120031 . hereafter BB03) find a notable difference between 
the luminosity functions of the core and halo of M67, with the halo having significantly 
more faint star contributions. This is found as well in several other studies, such as 



Montgomery et al. 



(119931 ) and 



Fan et al. 



(1996) 



2006 



hereafter 



BB03 make use of the Two Micron All Sky Survey (jSkrutskie et al. 
2MASS) to probe M67 over a large and continuous area. M67 has a large enough spatial 
projection on the sky that it could not be imaged by most large telescopes in a single 
exposure. Unlike dedicated "pencil-beam" surveys which could be carried out over 
individual p ortions of the clus ter, large scale surveys such as 2MASS and the Sloan Digital 



Sky Survey (I York et al. 



2000l . hereafter SDSS) can uniformly image regions much larger 



than a single cluster. This makes them ideal resources for studies of cluster structure. 

In this paper we use SDSS data to probe the spatial and stellar structure of M67 as 
it is seen today. Our selection of SDSS photometry is described in §2. The matched filter 
technique which separates M67 from the surrounding field population is outlined in §3, and 
our results are presented in §4. In §5 we offer a discussion of the results and a calculation 
of the orbit of M67, as well as a summary to put our findings in context with previous work 
on M67. 
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2. SDSS Photometry 

The SDSS is a major survey consisting of both five-band (ugriz) photometric and 
optical spectroscopic data. The photometric data are collected in adjoining stripes over 



a quarter of the sky, covering the nort h galactic cap. It 



ranging from near by asteroid families 



Parker et al. 



las been used to study objects 



(llnada et al. 



20081 ) to incredibly distant quasars 



20091 ). The SPS S provides nearly simultaneous i maging in five photometric 



bands (ugr iz), down to r~23 



calibrated ( ITucker et al. 



(Adelm an- McCarthy et al. 



2006 



Davenport et al. 



20081 ). These data are carefully 



20071 ). achieving excellent photometric 



errors averaging 1% in u,g,r,i and 2% in z. The large spatial coverage combined with 
accurate photometry has been useful for studying low-density str uctures such as tidal 



Grillmair fc Dionatos 



stream s in the galactic halo and globular cluster disruptions (e.g. 
2006bl ). The SDSS thus provides an ideal source of information about several nearby open 
clusters, which are often distributed over large portions of the sky. 



Our data came from the sixth SDSS data release ( Adelman-McCarthv et al 



2008 



hereafter DR6) which is publicly available online via a SQL databasqj. The 
reader is referred to DR6 and references therein for a highly detailed description of 
the survey. Our query returned sources from the STARS view of the PhotoO biAll 



20011 ) detects 



table, which includes all objects which the SDSS pipeline (ILupton et al. 
as point sources (as opposed to extended sources or cosmic rays for example) and 
are not quasars. To ensure a clean sample of only stars in our analysis we set the 
following photometric flags: SATURATED = 0, BRIGHT = 0, OK_SCANLINE > 0, 
OK_STRIPE > 0, EDGE = 0, LOCAL_EDGE = 0, PRIMARY > 0, PSF_FLUX_INTERP 
= 0, INTERP_CENTER = 0, BAD_COUNTS_ERROR = 0, PEAKCENTER = 0, 
NOTCHECKED = 0. Some of these flags are redundant to those automatically applied 



] http : / /cas j obs . sdss . org] 
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when querying from the STARS view. The stripes of the survey run across great-circles 
along galactic latitude. We found flags such as OK_SCANLINE, OK_STRIPE, and 
PRIMARY particularly necessary to remove over-densities of stars found within the 
overlaps of the stripes. The SDSS provides a description of suggested use of flags online at 
http : //www, sdss . org/dr6/products/catalogs/f lags ,~html| We have chosen to use the 
PSF magnitudes as they are considered the most unbiased measurement for point sources 
by the SDSsE 

We queried a 20° x 20° box centered around M67 (a=132.825, £=+11.8, J2000) which 
returned 2,111,236 point sources with 15 < g < 23 and 14 < r < 22.5. These magnitude 
cuts were made to remove major photometric scatter which grows rapidly near the faint 
limit, and to avoid bright stars which may have spurious photometry. The entire sample is 
shown in Figure [TJ which shows the density of point sources found by SDSS as filled contours 
(where darker is increasing density 



SDSS from the 



Schlegel et al. 



Estimated photometric reddening is computed by the 
(ll998f ) dust maps. We carried out the analysis in the following 
sections both with raw magnitudes and reddening corrected magnitudes. These reddening 
adjustments are critical for extragalactic or distant cluster science as they correct for the 
extinction of light from the interstellar medium. Because M67 is quite close, extinction is 
not suspected to be a major contribution to the observed photometry. Indeed, no significant 
differences arose in our results when using de-reddened magnitudes. The final analysis was 
completed with the de-reddened data. 



2 http : //www . sdss . org/DR6/ algorithms/photometry . html 
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The Matched Filter 



The matched filter technique, als o known as optimal contrast filtering, has been used 

hereafter R02) used a matched 



in galactic astronomy for some time. 



Rockosi et al. 



(120021 



filter on early SDSS photometry around the globular cluster Palomar 5. They were able to 
detect large symmetric tails of stars leading and trail ing the cluster in its o r bit. Fo llow-up 



matched filter work by 



Odenkirchen et al. 



( 120031 ) and iGrillmair fc Dionatosl (j2006al ) showed 



these tails to extend at least 22°, a limit imposed by the edge of the SDSS footprint. We 
used the DR6 data around Palomar 5 as a benchmark for our algorithm, and were able to 
recover the ~22° tails. A similar algorith m has been employed to de t ect low-density stream s 



of star s with no detectable nucleus (e.g. 



Grillmair fc Johnson 



2006 



Grillmair fc Dionatos 



2006b|). R02 provide a good description of the matched filter process, and the reader is 
encouraged to seek further discussion therein. We will however briefly describe our method 
and its application to the open cluster M67 as well as our artificial clusters below. 



3.1. Analysis of M67 from SDSS 

Traditionally in the absence of more detailed information, such as proper motions, 
radial velocities, or spectra, cluster membership can be crudely estimated based on the 
positions of stars in a color vs magnitude diagram (CMD). To maximize the contrast of the 
cluster against the surrounding field stars we must characterize the probability functions 
from both the field and cluster in color-magnitude space. These two CMD probability 
functions are then used to determine the probability of membership for every star. 

The M67 CMD was determined by selecting the stars within 0.25° of the cluster center. 
Figure [2] shows the M67 CMDs for g — i vs g and r — z vs r in our SDSS data. The main 
sequence is clearly visible in both panels, and a faint equal-mass binary sequence can be 
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seen above it. We chose to use the colors g — i and r — z because they better separate 
the M67 main sequence from the significant field contamination. The galactic field stars 
were crudely filtered out by rejecting any stars not within a few photometric a of a spline 
hand-fit to the main sequence and equal mass binary sequence. A 2-D gaussian smoothing 
algorithm was then applied to these rough cluster CMDs, creating a continuous CMD 
distribution for M67. 

The number of cluster stars in a given solid angle is described by the equation 
n ci — otfch where a is the number of stars in a region, and f c i(color,mag) is a normalized 
probability function for the cluster in the color-magnitude plane. By normalizing the 
continuous CMD distributions described above, we created the fa functions for M67. This 
was done for both (g — i, g) and (r — z, r) independently, and are shown in Figure El The 
non-uniform distribution shown in Figure [3] is a result of sampling f c i from the stars in the 
core of M67 which are not evenly distributed along the main sequence. R02 investigated 
whether the mass function of the cluster core would detract from the matched filter's ability 
to detect tidal debris that would generally have a different mass function. Their conclusion 
is that the method is robust against such biases, but we do bear this in mind in our analysis 
later on. 

The same task must be carried out for the field to describe the background filter. 
It is clear from Figure [1] that the field does not have a uniform density of sources in our 
sample. Instead we see the density of stars rise at lower galactic latitudes, near the galactic 
plane. To investigate the effect the changes in the field star population would have on our 
analysis, we sampled the background population in several regions. Figure H] shows the 
background (g — i,g) CMD in four quadrants of the field. Other than the number of stars 
in each quadrant, the differences between these four CMD samples is minor, and we believe 
the changes in the CMD will not greatly affect our analysis, as was also concluded in R02. 
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Thus we make the assumption that the background has a uniform CMD structure across 
the field studie d, and that the stella r dens ity may be described by a low-order surface fit, 
as described in iGrillmair fc Johnson! ( 120061 ) . 



The background filter is referred to as nb g (color,mag), which, given any area on the 
sky provides the number of background stars in each color vs. magnitude bin within that 
solid angle. It is valid to use rib g and not fb g here because the background is a continuous 
and relatively smooth population, and the number of stars in any area can be estimated. 
By using all of the stars around M67 with a radius between 1.5° and 8° (1,011,776) our 
background population is very well sampled. This annulus was chosen to avoid both the 
cluster and any potential tidal features, and the large hole in the data set in the south-east 
corner. Applying a boxcar smoothing algorithm and dividing by the area used, we create 
the two ribg filters shown in Figure El 

Because the field population does not have a uniform distance or age, the familiar 
features of an ideal CMD (e.g. main sequence, turn off, red giant branch, etc) become 



severely blurred, 
diagram (e.g. see 



he resulting color vs. apparent magnitude diagram is known as a Hess 



Alcock et al 



20001 ) . In the Hess diagrams in Figures H] and [5j several 
distinct features are visible in the CMD contours. The peak at (g — i,g) = (0.6, 17) is from 
the galactic disk (the thick disk according to R02). This feature is strongest in Figure HH, 
which has the lowest galactic latitude and thus should contain the most disk stars. The 
smooth feature at g — i ~ 0.5 is the main-sequence turnoff for the halo population, which is 
naturally most prominent in Figure Hh, at the highest galactic latitude of our sample. The 
large red buildup of stars at g — i ~ 2.5 is the local field K and M-dwarf population, most 
of which are likely part of the disk. This is also most prominent in Figure HH. We have 
enforced magnitude cuts at both the bright and faint limits of our SDSS sample, and thus 
the f c i and n& g distributions have steep declines near the limiting magnitudes. 
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To produce the optimal contrast, the normalized cluster filter is divided by the scaled 
background. This yields h — which is shown in Figure El This produces a CMD 
matched filter that describes the stars which stand out from the background the strongest. 
For any given star, its color and magnitude will yield a value of h related to its probability 
of membership in M67. Naturally this probability function will promote regions of the 
CMD where there are many M67 members but few field stars, and conversely punish stars 
with CMD positions having large field populations. 

The spatial distribution of the cluster is then mapped by summing the values of h 
inside small (0.1°, 0.1°) spatial bins. We complete this entire proces s independently for bo th 



the (q 



i,g) and (r — z,r) CMDs, and co-add the results, as per 



Grillmair fc Dionatos 



(|2006aJ). While the probability distribution h produces the maximum likelihood indicator 
for the stars, corrections for the field response to the h filter must be done. To solve for the 
actual number of cluster stars found in each spatial bin, a, the following formula from R02 



is used: 

fcl 



a 



f 2 

Jcl 



jili _ / f d d(color,mag) \ I \ l -^-d(color,mag) J> , (1) 
where — = h is summed for all stars in a spatial bin as mentioned above, the integral of 

f 2 

f c i accounts for the background response in a solid angle d£l, and the integral of is the 
signal-to-noise response of the filtering. This equation is the basis for our analysis with 
both the SDSS and 2MASS data sets, as well as the Monte Carlo simulations described in 
the following section. 



3.2. Analysis of Model Clusters 



In order to test the robustness of the matched filter method, we ran our algorithm 
using artificial clusters of known stellar composition and structure. By characterizing the 
efficiency and reliability of the algorithm we were able to examine issues of biases and 
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errors for our results in §4. To accurately measure the efficiency and determine the greatest 
sources of error, we inserted artificial clusters at random locations within our field. These 
tests were run 1000 times on our SDSS data set in order to determine the reliability and 
significance of low-contrast features seen in §4. 

The simulated clusters were created with a large range of core densities. The number 
of members was randomly chosen between 450 < Nstars < 1350, while the core radius was 
randomly selected between 0'.9 < R CO re < 22'. This wide range was used to fully explore the 
densities and spatial compositions which open clusters might have. 

The spatial composition of each cluster was created by randomly generating radial 
positions for every member, each being drawn from a gaussian profile with a standard 
deviation radius of R core . The angular positions were then chosen randomly between 0° 
and 360°, creating a roughly uniform circular distribution of stars. No tidal elongation 
was explicitly included, however the gaussian distribution of stars did at times produce 
asymmetric halos at large radii consisting of up to a few dozen stars. 

The griz magnitudes were formed using the handmade splines fit to M67 in §3.1. A 
luminosity function for the g band was used which had the form g = (g®' 77 ) x 8 + 15, where 
is an array of Nstars elements filled with random values between and 1. This function 
was chosen to provide an increasing number of stars with decreasing luminosity, but to 
have a deficiency as compared to the field star population of low-mass stars, as observed in 
old open clusters. We also found that it approximately reproduces the observed luminosity 
function of M67. An equal mass binary sequence was created by increasing the magnitudes 
of a subset of stars by 0.75 mag in every band. The fraction of equal mass binaries was 
chosen at random between 10% and 60%. 

Appropriate photometric scatter was created by sampling errors for the actual SDSS 
data. The mean error and the standard deviation of the errors were calculated for griz 
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bands in bins of 0.1 magnitudes. These errors were then mapped to the artificial stars 
by their corresponding griz magnitudes, choosing for each star a random number within 
the gaussian error envelope. Adding these errors to the artificial photometry created a 
very realistic scatter about the M67 splines which increased with increasing photometric 
magnitude. 

For every model run, the simulated data was placed randomly in the field around 
M67, but was required to be greater than 1.5° radially away from M67 to avoid cross- 
contamination, and less than 8° away to stay well within the bounds of the SDSS data 
sample. The full matched filter analysis code was then run, sampling f c i from a 0.25° radius 
around the model cluster and smoothed as before with the M67 f c \. The background n\,g was 
determined from an annulus around this with a radius between 1.5° and 6°, encompassing 
~500,000 stars on average, but taking care to avoid M67. 

For every run of the simulation we recorded the simulated cluster position, surrounding 
background density, the estimated number of cluster members recovered by our algorithm, 
the radial surface density profile for each run, as well as every plot used in the analysis of 
M67. In Figures [7| and [8] we show example results for a random subset of our models. The 
resulting statistical analysis of the model runs show that we are able to reliably recover 
both the number of members and their locations. 

Figure [9] shows the mean recovered radial surface density profile for our 1000 model 
runs, marked by stars. The squares show the mean radial surface density profile of the 
1000 input models themselves. We are on average over estimating the surface density in 
the central bins, but are recovering the core radius well. Triangles mark the mean residual 
between the input models and the recovered profile. The error bars are one standard 
deviation in each residual bin. The small overestimation of the radial profile in the nucleus 
is attributable to the spatial smoothing kernel which will indicate more stars in the central 
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radial bin than are actually observed. 

In every run of our models the radial surface density profile is fit with a least-squares 
gaussian profile. From the 1000 model profiles we found on average our algorithm 
underestimated the core radius by 11%, and after correcting this bias we found a standard 
deviation of 0'60 for R c . The number of stars detected by our algorithm was larger than the 
input models by, on average, 17%, and we computed a standard deviation on the number 
of stars recovered to be 67 stars. We also determined that the input binary fraction was 
1.27 times more than our recovered binary fraction, and found a standard deviation of 6% 
for the resulting binary fraction. None of the recorded statistics changed as functions of 
the background density, suggesting that our matched filter algorithm consistently was able 
to remove the field contamination. We are attempting to improve these biases for future 
implementation of our algorithm. 

Since each model run was required to be at least 1.5° away from M67, in some cases 
the artificial cluster can be close enough to M67 for it to be present on the radial profile 
plot at it's furthest extent. This can be seen in Figure [Sh^ for instance. These examples 
provide a useful benchmark for the robustness of the matched filter. Despite low-order 
changes in the local background population surrounding each run, and the variations in the 
model fd compositions and model core densities, M67 is always nearly perfectly recovered. 
The elongated features of M67 seen in these examples match those found in the analysis of 
the following sections. 

Visual inspection of a subset of the models, as demonstrated in Figure El showed no 
examples of spurious extra-tidal features which deviated from the gaussian distribution 
caused by fluctuations in the background population. Any elongated features appeared to 
be caused by the input model. Further, our algorithm itself, while slightly overestimating 
the number of stars in the clusers' core, does not produce random cluster halo signatures. 



- 15 - 



We therefore conclude that the matched filter method is effective in characterizing and 
removing the field star population and providing an accurate representation of the cluster 
within the range of our sample. Any extra-tidal structures or elongations found in the 
vicinity of M67 are therefore considered to be real and intrinsic to the cluster. 

4. Results 

4.1. M67 Properties 

By analyzing SDSS data for M67 using the technique outlined in §3.1, we were able to 
map the cluster distribution to a larger radial extent than has been done before. Figure [TU1 
shows our spatial map of the open cluster M67. A small (0.1°) boxcar smoothing kernel 
has been applied to the summed a data, and a low-order surface fit to the background was 
subtracted to remove any residual large-scale variation of the field density. The core of 
M67 is strongly detected, and shows a circular distribution. The over-plotted circle, 25' in 
radius, denotes the furthest radius which BB03 detected the cluster against the surrounding 
field population using 2MASS data. There is however a significant low density asymmetric 
halo of stars well outside the core. The contours from light to dark are increasing levels of 
density. They are defined by the equation levelj = MED(a) + (a Q x 2 J ) where MED(a) 
is the median value of a over all spatial bins, a a is the standard deviation of a, and 
j = {1,2,3,4,5,6}. To ensure that a was well characterized, we created a histogram of 
a from every 0.1° x 0.1° spatial bin in our SDSS sample. This is shown in Figure [UJ 
and appears gaussian in shape and is centered around a = 0. Since our contours begin 
at two standard deviations above the peak shown in Figure [TTl and noting the reliability 
determined in §3.2, we argue that these features, discussed in §4.2, represent real structure 
around M67 which has previously gone unseen. 
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The radial surface density profile of M67 is shown in Figure [12j The strength of 
detection in each radial bin is the sum of a values in that bin, divided by the area that the 
bin encompasses. This is directly analogous to the spatial map in Figure [TD1 and the units 
of (a arcmin~ 2 ) are equivalent to (stars arcmin~ 2 ). The extended halo seen in Figure [TUl 
is visible from a radius of ~ 25' to ~ 50' where it falls to background levels. 



19621 ) to our data, 



We fit this surface density profile using a King-like profile ( iKingl 
which has the form n(r) = n^kgd + ^ , where R c is the core radius, n bkgd the 

background surface density (not to be confused with the background filter n& s ), and no the 
central peak surface density. This surface density profile equation was used perviously by 
BB03 to model M 67. This fit is shown in Figure [12] as the solid lines. Our King model fit 
provides R c = 8'24 after correcting for biases from §3.2, somewhat larger than the previous 
2MASS based measurement of R c = 4f86 by BB03. From our testing with the model 
cluster, we estimate the error to be 0'.60. We believe the difference between our radial 
profile and that of BB03 is due primarily to the lower limiting mass our study probes. We 
do give a direct comparison to BB03 by using the matched filter algorithm on the M67 
2MASS data and provide more discussion in §4.3. After fitting the first-order King model 
used in BB03, we attempted to determine the tidal radius for M67, R f . Usi ng the equation 
n{r) = n bkgd + n Q {[l + (r / 'R c ) 2 ]- 1/2 - [1 + (R t / R C ) 2 Y 1/2 } 2 from|King| \^M) and the values 
we determined for the first-order fit above, we found the estimated tidal radius for M67 to 
be R t = 64' ± 21'. Errors here are determined from standard deviation of the RMS scatter 
on our fit. This corresponds to a radius of 16.8 pc at the distance of M 67, much larger than 



the determination using bright stars in M67 by iPiskunov et al. 



( 120071 ) . By summing the 



number of stars in each radial bin in Figure [121 we estimate the total number of visible M67 
members to be 1385 ± 67, where our error is again adopted from our models in §3.2. 
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4.2. Composition of the Halo 

Having found evidence for a large halo around M67, we investigated the differences 
between the halo and core populations. Variations in the types of stars in the inner and 
outer regions of the clusters might suggest an origin for the elongated halo. In dynamically 
un-evolved systems, the initial stellar luminosity function rises towards lower masses. BB03 
use this property to investigate the evidence for mass segregation in M67. They created 
a set of 2MASS J-band luminosity functions in three regions (core, inner-halo, and outer 
halo). A higher fraction of low- mass and faint J-band stars were seen in the outer halo, and 
thus mass-segregation is implicated. 

Binary star systems are composed of two stars which are often too close to be visually 
resolved. This single point source is observed to be brighter than a single-star of the same 
apparent color. These systems also have higher mass than single-star systems by definition. 
The effects of mass segregation discussed above are therefore applicable to binary stars 
which act as more massive single stars. Previous cluster studies have observed a stronger 
concentration o f binary stars in the cores than in the outer envelopes of older clusters (e.g. 



Fan et al. 



1996|). 



The binary fraction and luminosity function for each region of our data can be probed 
in much the same way as the spatial structure. Rather than summing a, the number of 
potential member stars, in each spatial bin, we summed a in three concentric regions for 
every magnitude bin (for the luminosity function) or across the main sequence (for the 
binary fraction). These three regions were the core, halo, and a large background annulus 
centered around M67. 

Calculating the binary fraction across the entire magnitude range requi res properly 



collap sing the CMD into the color plane. We employed a technique used by 



Clark et al. 



(120041 ) to collapse the main sequence of the globular cluster Palomar 13. Since a spline 
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has already been fit by hand along the M67 main sequence in §3, we subtracted the 
main sequence spline color from every star's color, thus centering the main sequence on 
AColor = 0. Because the equal-mass binary sequence is 0.75 mag brighter than the main 
sequence, it's AColorbin position was easily tracked. We then divided AColor by the 
AColorbin of the binary sequence, creating a reduced color R = AColor / AColorbin- This 
placed the main sequence at R = and the equal mass binary main sequence at R = 1. To 
calculate the binary fraction we summed all of the a values in each reduced color bin, with 
a range of —1 < R < 2. The luminosity function is found by summing all of the a values in 
every magnitude bin, again cutting out stars outside the reduced color range — 1 < R < 2. 

In Figs. [131 El and [15] we investigate the M67 core, halo, and background populations 
respectively. The central panel shows the reduced CMD, where the darker pxels show 
increasing a sums in each CMD bin. The right panel shows the sum of a values along 
the magnitude axis, while the bottom shows the same along the reduced color axis. Since 
the equal-mass binary main-sequence was created by translating the main-sequence spline 
0.75 magnitudes brighter, a discontinuity in the calculation of AColorbin arises at the end 
points. Thus we have reduced our magnitude range by 0.75 on both the bright and faint 
end to avoid these biases. 

The core population in Figure [[3] (radius < 0.8°) shows a clear double gaussian peak 
along the reduced color R axis, and a deficiency of faint (low-mass) members along the 
magnitude axis. The halo in Figure HH (0.8 < radius < 1.2°) contains a more flattened 
luminosity function, and a diminished secondary peak in the binary fraction. These two 
figures do not however properly account for the contribution the field population would 
have, especially in the halo which samples a larger area and has a lower expected density 
of cluster members. Thus in Figure [15] we made the same measurements for a much larger 
sample of stars (1.5 < radius < 6°) to estimate the field contamination in our binary fraction 
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and luminosity function. This shows the response to the h filter that the background has. 

The background contribution was normalized and then scaled by the area in the 
core and halo regions respectively. This bias contribution was then subtracted from the 
M67 core and halo samples. Figure [16] shows the background subtracted binary fraction 
estimations for the core and halo of M67. A gaussian profile was fit to the primary peak for 
R < 0.25. The residual in the core was nicely fit by a second gaussian profile, centered at 
R « 1. Integrating the two gaussian profiles, we measure a fraction of bin ary point sources 
to be 21% in the co re of M67. This is c ompa rable to the 



Fan et al. 



16% and 22% from 



Montgomery et al. 



(119961 ) determination of 



(ll993f ). Correcting for the artificial cluster tests in 



§3.2 above, our estimation of the binary fraction is increased to 26%. This is necessarily a 
lower limit on the fraction of binary sources as many stars may be binary systems without 
detectable flux excess. No significant gaussian binary star residual was found in the halo 
population, although the signal is noisier. Both the core and halo samples go to a = for 
R < and R > 1. The luminosity function is given the same treatment, and the background 
corrected functions are presented in Fig. [TTJ The halo shows a more flattened low-mass 
contribution. The core however contains a significant lack of low-mass members, even after 
accounting for the background field contribution. 

The dichotomy between the core and halo populations is consistent with qualitative 
ideas about mass segregation. As discussed in §1, an old cluster such as M67 is expected 
to be losing low-mass members to the outer regions of the cluster, while the core becomes 
increasingly more concentrated with higher mass stellar systems. This simplification of 
course assumes a single static population of stars and also ignores the effects of classical 
stellar evolution, whereby higher mass stars die before lower mass, and skew the observed 
present day mass function. Assuming an age for M67 of ~4 Gyr, it would be expected that 
late A or early F stars would be leaving the main sequence, and thus would be brighter 
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than our limiting magnitudes with SDSS. We do not therefore expect stell ar evoluti o n to 
contribute significantly to the observed luminosity functions of our sample 



Anetal 



(120081 ) 



investigated the effects of crowding and stellar density on the photometric completeness 
of the SDSS pipeline, including the M67 region. They found that the automated pipeline 
recovered a comparable number of photometric sources as compared to manual reductions of 
the imaging for the open clusters M67, NGC 2420, and NGC 6791. Thus we conclude that, 
while manual reductions of wide-field imaging such as in our sample would be preferable, 
the SDSS pipeline is able to reliably probe the regions around high galactic latitude open 
clusters. 



4.3. Comparison to 2MASS 

To provide a complete comparison of our matched filtering method with previous 
results for M67, we ran our algorithm on J, H, and K s band photometry from 2MASS. Using 
the 2MASS interface available on the internet entitled GatoJ^, we retrieved a 10° x 10° 
box surrounding M67. This produced 264,354 point sources. The cluster CMD probability 
distributions f c i(J — H, J) and f c i(J — K s , J) were drawn from all stars within a radius of 
0.25° around M67, and the background CMD from everything with a radius greater than 
1.5°. Our IDL code for the SDSS data set was modified to use the 2MASS JHK S data, 
and the analysis was carried out to measure the spatial distribution of M67. We did not 
measure the binary fractions or luminosity functions for the halo and core populations with 
the 2MASS data. 

M67 as seen by 2MASS contains a different range in spectral types than is observed by 
the SDSS, reaching significantly higher masses. This is shown Figure 2 of BB03, with the 



3 http:/ /irsa.ipac. caltech.edu 
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M67 CMD reaching a spectral type KO at the faint limit, which we adopted to be J = 17. 
In addition to the upper main sequence, the 2MASS CMD also contains the turn-off and red 
giant branch. The average stellar density in our 2MASS sample was 0.73 stars arcmin~ 2 , 
whereas the average in our SDSS sample was 1.4 stars arcmin -2 over the same spatial area. 

Figures [TS] and HH summarize our 2MASS results for M67. These are analogous to 
Figures [10] and [12] respectively. An asymmetric distributions of stars are seen in Figure 
[18] which extend beyond the detection of the BB03 study. The elongation seen in Figure 
[TO] is not well reproduced by the 2MASS data, however the mass ranges are considerably 
different. 

Two significant differences between the SDSS and 2MASS results are apparent in 
Figure HH] the 2MASS data yields a significantly higher surface density in the core, and the 
core radius fit with a King profile is much s maller (R r = 4? 12 ) . This core radius corresponds 



very well with determinations by BB03 and 



Fan et al. 



(119961 ). Because the 2 MASS sample 



contains higher mass stars than the SDSS sample, the discrepancy with §4.1 in core radius 
is expected. The higher mass stars in the 2MASS data will be more centrally concentrated 
in the core. These are too bright to be observed in the SDSS data and thus the SDSS will 
probe a mass range with a larger core radius. 

Mass segregation differentiates all of the stars in a cluster according to their mass. The 
SDSS and 2MASS core radius comparison provides two estimates with a large separation 
of mass. To test that mass segregation was still evident in the 2MASS data alone we split 
our 2MASS sample into a bright bin and a faint bin. The separation was chosen to be 
J = 12.5 which roughly corresponds to the main sequence turn-off. The full matched filter 
was rerun for each 2MASS subset. The bright and faint samples yielded core radii of 3'70 
and 4'.Q7 respectively, indicating that mass segregation among the high mass sample is 
clearly evident, as found in §4.2 above. 



-22 - 



The asymmetric halo arou nd M67 shown in Fig. [TO] appears elongated roughly in-line 
with the proper-motion vector. 



Fanetal 



( 119961 ) have previously suggested an elongation 



of the core in M67 which roughly matches what is seen in Fig. [TS] but d oes not resemble 
what is seen in Fig. [TUJ A str ong correlation c annot be made between the 



Fan et al. 



results and our work because 



Fan et al. 



(1996) 



119961 ) considered asymmetries through 8 angular 
bins covering 45° each. Their optical data was also limited to a photometric depth of V = 21 
and 0.3 magnitude errors makin g comparison to our SDSS data tenuous. Qualitatively the 



core radius and tidal features of 



Fan et al. 



(|1996I ) match our 2MASS analysis of M67. Our 
matched filtering has been able to reproduce known structure for M67 using 2MASS, and 
revealed a slightly increased radius of detection for the halo compared to BB03 due to our 
more complete subtraction of the contaminating field population. 



5. Discussion and Summary 

We have used a matched filter algorithm on SDSS photometry for identifying probable 
members of the open cluster M67. This study has revealed a core radius R cor e = 8f24±0f60, 
and an asymmetric distribution of stars outside the previously known tidal extent. This 
asymmetric halo extends nearly 60' from the cluster center. The total number of stars 
measured within our SDSS data for M67 is 1385 ± 67. 

We have also tested the reliability of our matched filter method using 1000 artificial 
cluster simulations. In our models we tended to overestimate the number of members by 
17% on average, and underestimate the core radius by 11%. Overall the matched filter 
technique appears to be a robust method for identifying low-density stellar populations 
against a significant background. 

Our measurement of the binary fraction in the core of M67 is in good agreement 
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with many previous studies (e.g. Montgomery et al.lll993t iFan et al.lll996l ). Correcting for 
systematic biases measured with our models, we estimate a lower-limit on the fraction 
of binary sources detected to be 26 ± 6%. Since each binary source is composed of two 
stars, this corresponds to at least 45% of all stars in our M67 core sample belonging to 
binary systems. This represents a lower limit as many of the "single" stars are possibly 
bound to very low mass stars which would place them blue-ward of the equal mass binary 
main sequence. We als o do not correct for stars which may be obscured within the cluster. 
Bica &: Bonattol (120051 ) measured a binary fraction for the core of M67 to be 39 ± 16%, 
which agrees with our determination within the respective errors. The increased central 
concentration of binary stars in the core, along with the greater fraction of low mass stars 
in the halo, fulfills the expectation by several studies that M67 has indeed underg one 



Hurley et al. 



20051 ) are to 



significant mass segregation. If N-body simulations for M67 (e.g. 
be believed, this cluster is near the end of its life as a gravitationally bound and distinct 
stellar population. 

BB03 have recently carried out a very similar examination on M67 using 2MASS. 
The relatively bright limiting magnitude in the infrared meant that BB03 sampled higher 
mass stars than the SP SS. Our SDSS DR6 sample for M67 reaches a spectral type of M4 



(iBochanski et al. 



20071 ). but is limited at the bright end to mid K type stars. It is expected 



then that our SDSS sample would show signs of mass segregation to further radii, as we 
probe significantly lower masses. Since the elongated halo is not reproduced with our 
analysis of the 2MASS data, we suggest it must be primarily made up of the lowest mass 
members in M67. 



Bergond et al. 



(120011 ) have suggested that a general alignment between the proper 
motion and the elongation of an open cluster is not uncommon. They also demonstrate 
that the major axis in the core of a cluster may preferentially point toward the galactic 
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center. The second highest contour in Figure [10] in the core of M67 is slightly elongated 
along the direction of the galactic cent er. Our newly discove red tidal tails also emanate 



almost perpendicular to this as seen in 



Bergond et al. 



(1200 ll ). and our tidal features seem 



to fit all expectations of an open cluster being pulled apart. 



Carraro fc Chiosil (119941 ) project the orbit of M67 back in time to approximately its 
birth. Their calculations show that M67 has made approximately 17 passes through the 
galactic plane [Z — 0), and that the current position of M67 is near the vertical apex of its 
orbit (Z = 0.41Kpc). It is generally thought that the large distance above the galactic plane 
at which M67 spends much of the time has helped prevent tidal effects from completely 
destroying the cluster, as is the case with NGC 6791. 

Usin g modern determinations of the proper motions, radial velocity, and distance to 



M67 from 



Frinchaboy fc Maiewskil (|2008r). we ca 



followed the same procedure as 
galactic potential, and used the 



Carraro fc Chiosi ( 



culate d the cluster's full galactic orbit. We 



1994]), which employed a three component 



Johnson fc Soderblom 



space velocities. Figure 12171 shows the meridional projection for the orbit o: 
5x 10 8 years into the past and future. This is in excellent agreement with 



(119871 ) method of computing the 

M67, calculating 



Carraro &: Chiosi 



(11994] ) and earlier estimations. 



Typically long streams of tidal debris, such as those trailing /leading Pal 5, are assumed 



to lie along the orbital path of the cluster dMontuori et al 



has re cently come under scrutiny for Pal 5 (jOdenkirchen et al. 



2007) . How ev er this assumpt ion 



2009|). 



Bergond et al. 



(120011 ) also show examples of tidal debris around open clusters which is not grouped into 
dynamically cold streams. In order to determine the projected direction tidal debris tails 
would be expected to lie along for M67, we computed star paths in a reference frame moving 
with the average orbital velocity of the cluster. Simply put, stars that are leading or lagging 
the main body of the cluster in its orbit around the galaxy do not diffuse to completely 
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ring the galaxy because they mostly share the orbital motion of the cluster. Rather they 
will occupy a more modest range of azimuthal angles in the galaxy, with smaller ranges 
corresponding to smaller orbital eccentricity. M67 moves both in the radial and vertical 
directions (as shown in Fig. 120]) . For an open cluster orbiting with vertical excursions from 
the disk, the orbit will not be closed, and the cluster (and dynamically cold stars that have 
escaped the cluster) will eventually sample different parts of a boxy volume moving with 
the average orbital velocity of the cluster. To visualize this three-dimensional dynamical 
structure, we subtracted the average angular velocity of the cluster multiplied by the time 
since the present day from the azimuth angle. The projection of this motion onto our line of 
sight towards M67 looks nearly straight, as shown in Fig. [10] Any kinks or sharp changes 
in the projected co-moving deviation are well outside our SDSS field. 

Since the most elongated feature of M67's halo seen in Fig. [TD] is not aligned with 
the projected co-moving deviation, or the projected orbital path of M67, we believe the 
asymmetric tidal feature's origins seen with SDSS cannot be singularly attributed to the 
galactic tidal field. However a separate smoking gun, such as a nearby molecular cloud or 
cluster, is not readily apparent to us. It is more probable that with such small numbers of 
stars found outside the core, tidal shocking from passing through the galactic disk produces 
these types of weakly elongated features in open clusters, rather than the dramatic tails as 
seen in Pal 5. This can only be verified, however, using a detailed survey of the kinematics 
of the cluster members on a large spatial scale to compare the core and halo dynamics, 
along with future N-body modeling of cluster disruption. 



The optimal contrast filtering we have employed here will be of 



the n ext generation of wide field surveys. Programs such as LSST ( llvezic et al. 



great use in 



2008|) and PanSTARRS flKaiser et al 



20021 ) will map the sky with never before seen 



levels of photometric precision and depth, with a spatial coverage far exceeding even 
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SDSS. The most recent SDSS public release ( lAbazajian fc Sloan Digital Sky Survey 



2008 



DR7) includes the Sloan Extension for Galactic Un derstanding; and Exploration 



( INewberg k, Sloan Digital Sky Survey Collaboration 



20031 ) data, which includes photometry 



in and around the galactic plane. This new data set includes well over a hundred open 
clusters, with a great range in ages, masses, and surely dynamical states. We anticipate the 
matched filtering technique will help detect many new open and globular cluster features, 
as well as continue to find other faint substructure in the galactic halo. 
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Fig. I.— The 20° x 20° SDSS DR6 field around M67. Darker contours indicate a higher 
density of stars. The stellar density ranges from roughly 3300 to 11100 per deg 2 , with the 
density near M67 at ~5000 per deg 2 . Contours are in increments of 2000 stars per deg 2 . M67 
is faintly visible at the center of the field, and increasing galactic coordinates are marked 
from its center. The swath of missing data starting at about (a = 123, S = 5) is due to 
a missing piece of one stripe of the survey. Roughly a dozen white spots are seen, which 
mark locations of bright stars or galaxies which were masked out of the point-source catalog, 
including a bright star near M67. 
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Fig. 2.— The (g - i,g) and (r - z,r) CMDs of all 1389 stars within a radius of 0.25° from 
the center of M67 in our data. 
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Fig. 3. — The normalized and smoothed CMD filters M67, fd(g — i,g) (left), and f c i(r — z,r) 
(right). The contours are steps of la above the mean level of the filter. 




0.0 1.0 2.0 3.0 0.0 1.0 2.0 3.0 



Fig. 4. — The background filters n\,g{g — i,g) of four 5° x 5° regions in our sample. The 
(RA,Dec) centers of each region are a: (138.5,17.5), b: (127.5,17.5), c: (138.5,5.5), d: 
(127.5,5.5). Region (d) greatly overlaps the missing data seen in Figure [lj but the qual- 
itative features seen in panel (d) above are sufficient for our needs. Contours in all four 
panels are in steps of \o a above the mean level of panel (a) to demonstrate the increasing 
stellar density at lower galactic latitudes. Specific features are discussed in the text. 
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Fig. 5. — The smoothed background CMD filters, n bg (g — i,g) (left), and n bg (r — z,r) (right). 
The sample is drawn from a region around M67 with 1° < radius < 8°. The contours are in 
increments of \a above the mean. Specific features are discussed in the text. 
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Fig. 6. — The matched CMD filters, h(g — i,g) (left), and h(r — z,r) (right). The contours 
are la increments above the mean of the h functions. 



Fig. 7. — The radial surface density profile for a random subset (4 of 1000) of our artificial 
cluster model tests. The solid line is a gaussian fit to the data, while the dashed line is a 
King-like profile fit. The dark stars are the recovered density at each bin, and the squares 
are the input model for this run. There is generally very good agreement with the input data 
to the shape and distribution of the measured results, although we tend to overestimate the 
number of stars in the central bins. 
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Fig. 8. — The spatial map of our model cluster for a random subset (4 of 1000). These 
correspond to the radial surface profiles in Figure Contours are in increasing a (cluster 
stars), as defined for Figure [TOl 
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Fig. 9. — The average radial surface density profile recovered from all 1000 of our models, 
shown as stars. The boxes are the average surface profile for the 1000 input models. Triangles 
represent the average residual between the input model and the recovered profile. Error bars 
are the standard deviations in each residual bin. 
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Fig. 10. — The spatial map of high probability cluster members for M67. Darker contours 
represent increasing exponential density of members, a. The contours are defined by the 
equation level j = MED(a) + (a a x 2 J ) as defined in the text. This produces the first contour 
at 2a. Galactic (I, b) coordinat e arrows and the corrected proper motion vector (PM) from 



Frinchaboy &: Majewskil (120081 ) have been added for reference. The arrow labeled g.c. points 



towards the galactic center. The dark circle denotes the furthest cluster membership in 
BB03, approximately 0.5° in radius. For scale, the bar measures the apparent size of lOpc 
at M67, a distance of 900pc. The dashed line shows the projection along the orbit which 
tidal debris would be expected to be scattered. 
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Fig. 11. — The resulting background a counts, in stars per 0.01 deg 2 . The entire sample 
has had a low-order polynomial surface fit subtracted which removes both large scale stellar 
density gradients (arising from known galactic structure components) and the bias from the 
response of the field stars to the matched filter. 
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Fig. 12. — The radial density distribution of stars around M67. The vertical axis is the sum 
of a values inside each radial bin. The solid line is a King-like profile fit which is described 
in the text. The dotted line is the fit to the data past 60'. Error bars are poisson \fN errors 
from the number of stars divided by the area of each radial bin. The insert is the same data 
shown in log-log space which is often used for cluster density profile studies. 
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Fig. 13.— The core of M67 (radius < 0.5°). A reduced color (R g _ h g) CMD is shown in the 
center panel, where the pixels represent the sum of the h values for every star in each CMD 
bin for this region. These h values are summed over magnitude space (right) showing the 
measured luminosity function, and reduced color space (bottom). The bottom panel can be 
used to estimate the number of stars in the main sequence (centered at R — 0) and equal 
mass binary main sequence (R — 1). 
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Fig. 14. — Same as Figure IT31 for the halo around M67 (0.5° < radius < 1.2°). 
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Fig. 15. — Same as Figure fT3l for the background around M67 (1.5° < radius < 8°). 
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Fig. 16. — The relative binary fraction in the core (left) and halo (right). The single-star 
main-sequence is centered at R = 0, while the equal-mass binary main-sequence is at R — 1. 
Note the lower recovery of binary stars in the halo. 




Fig. 17. — The relative luminosity functions in the core (left) and halo (right). The core 
is deficient of lower-mass stars compared to the halo due to ongoing energy equipartition 
within the cluster. 
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Fig. 18. — The spatial map of high probability cluster members for M67 using 2MASS. This 
figure has the same contours and labeling as Figure [TTJJ 
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Fig. 19. — The radial distribution of stars around M67 with 2MASS. This figure has the 
same fit and error information as Figure [121 
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Fig. 20. — The meridional plane orbit of M67, projected lGyr in the past and future. The 
present location of M67 is denoted by the dark circle. 



